Hybrid soft-mode and off-center Ti model of barium titanate 
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It has been recently established by NMR techniques that in the high temperature cubic phase 
of BaTiOs the Ti ions are not confined to the high symmetry cubic sites, but rather occupy one 
of the eight off-center positions along the [111] directions. The off-center Ti picture is in apparent 
contrast with most soft-mode type theoretical descriptions of this classical perovskite ferroelectric. 
Here we apply a mesoscopic model of BaTiOa, assuming that the symmetrized occupation operators 
for the Ti off-center sites are linearly coupled to the normal coordinates for lattice vibrations. On 
the time scale of Ti intersite jumps, most phonon modes are fast and thus merely contribute to an 
effective static Ti-Ti interaction. Close to the stability limit for the soft TO optic modes, however, 
the phonon time scale becomes comparable to the relaxation time for the Ti occupational states 
of Ti„ symmetry, and a hybrid vibrational-orientational soft mode appears. The frequency of the 
hybrid soft mode is calculated as a function of temperature and coupling strength, and its its role 
in the ferroelectric phase transition is discussed. 

PACS numbers: 77.84.Dy, 77.80.Bh, 63.70.-fh 



I. INTRODUCTION 

Barium titanate is a typical representative of per- 
ovskite ferroelectrics, which undergoes a cubic-to- 
tetragonal structural phase transition at the Curie tem- 
perature Tq — 403 K. The soft-mode nature of this tran- 
sition has been determined by neutron scattering as 
well as by hyper- Raman scattering however, some 
open questions concerning the role played by the Ti ions 
have so far remained unanswered. Chaves et al. pro- 
posed a thermodynamic model following the assumption 
of Comes et al. 0| that the Ti ion occupies one of the 
eight equivalent off-center sites along the [111] directions, 
which predicts a nonzero microscopic dipole moment of 
each unit cell and hence a transition of an order-disorder 
type. 

The off-center displacements of Ti ions in the high- 
temperature cubic phase of BaTiOa have recently been 
confirmed by NMR experiments Jsj, which revealed the 
order-disorder dynamics of Ti ions to be coexisting with 
the observable displacive features of the TO soft mode. 
This immediately raises the question about the appropri- 
ate theoretical model for BaTiOa and related systems, as 
the widely accepted soft-mode description implies a cen- 
tral position of the Ti ion in the high temperature cubic 
phase, leading to a phase transition of a purely displacive 
type. Also, first-principles calculations of the electronic 
structure p do not support the idea that the ferroelec- 
tric distortion is due to the Ti ion 'rattling' in the oxygen 
cage. Some general features of the two types of phase 
transition have been discussed by Aubry Q on the ba- 
sis of a linear coupled double-well model. The crossover 
between the order-disorder and displacive transition was 
investigated by means of molecular dynamics calculations 
by Stachiotti et al. [1] who used a two-dimensional shell 
model of oxide perovskites. A microscopic tretament of 
perovskite ferroelectrics has been carried out by Girsh- 
berg and Yacoby (GY) who describe the degrees 



of freedom associated with the off-center displacements 
in terms of Ising pseudospin operators and introduce a 
linear pseudospin-phonon coupling. They derive an ef- 
fective coupling between the off-center ions, which leads 
to a pseudospin ordering transition. The corresponding 
transition temperature is shifted from the instability tem- 
perature of the TO soft mode. 



In this paper, we make an attempt to develop a meso- 
scopic treatment which combines the order-disorder and 
displacive features of the phase transition. We introduce 
a simple mechanism at the mesoscopic level, which is 
closely related to the above GY model, i.e., we adopt 
the linear coupling between the symmetrized occupation 
operators for the off-center Ti sites and the normal coor- 
dinates for lattice displacements. Assuming that on the 
time scale of Ti intersite jumps most phonon modes are 
fast, we show that these modes essentially contribute to 
a renormalization of the static Ti-Ti interactions, which 
can be either ferro- or antiferro-distortive. Meanwhile, 
close to the Ti ordering transition the time scale of the 
soft TO mode becomes asymptotically comparable to the 
characteristic time for the relaxation of the Ti dipole 
moment. Thus a new vibrational-orientational hybrid 
mode appears, which retains some of its original soft- 
mode character, however, its frequency and width is de- 
termined by the the relaxational dynamics of the Ti sub- 
system. In addition to the hybrid mode, a purely relax- 
ational mode describing the intersite jumps of the Ti ions 
exists. It should be noted that the unperturbed TO soft 
mode describes the anharmonic vibrations of all the ions, 
and thus also include the oscillations of Ti ions around 
their average equilibrium positions. Apart from the time 
scales involved, the conditions for the formation of a hy- 
brid soft mode may also depend on the strength of the 
Ti-phonon coupling, which thus plays the role of a control 
parameter in the present model. 
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II. COUPLED TITANIUM-PHONON SYSTEM: 
STATICS 

Following Chaves et al. 0] we define the occupation 
probabilities nu = {1,0} for the off-center Ti sites in the 
i-th unit cell, where Z = 1, 2, • • • , 8 according to the con- 
vention specified in Ref. 3]. Obviously, '^ina = 1. Next, 
we introduce a set of symmetry adapted linear combina- 
tions of the nu variables, which transform according to 
the irreducible representations Aig, Aiu, Tin, and T2g of 
the cubic group. 
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Here we have omitted the cell index i. The variables YiT, 
where F = 1, 2, • • • ,8 labels the symmetries in the above 
order, satisfy the relation 

Y^r = 1 • (5) 

This implies, for example, that the Ti,j polar modes (3) 
can be effectively represented by three independent Ising- 
type variables 0- 

We will assume that there exists a direct coupling be- 
tween the symmetrized occupation probabilities or 'pseu- 
dospins' YiY due to long range dipolar interactions, which 
has the form 

^"i-^-lEE^r^^rV, , (6) 

In addition, we consider the interaction between thepseu- 
dospins and the phonon normal coordinates Q^p |9|, [lOj , 

T^mt = — j= E E ^'iP exp(ig • Ri) , (7) 

^ i qpT 

where q is the wave vector and p the branch index of lat- 
tice normal modes, /J^ the coupling constant, and Ri the 
lattice vector. Introducing the phonon momenta P^p and 
frequencies w^p, we can write down the phonon Hamilto- 
nian 

'^ph ^ 2 E i^%Q?pQ-?P + PqpP-qp) ■ (8) 

qp 



It has long been established that in BaTiOs and other 
pcrovskite ferroelectrics il^] a soft TO phonon mode ex- 
ists. Here we assume that the soft TO mode has a wave 
vector qo and symmetry F = Ti^^s, where s is the corre- 
sponding branch index. The frequency of this soft mode 
vanishes at the stability limit To in accordance with the 
Cochran relation 

u;l,=a{T-To) . (9) 

In the following we will limit ourselves to the case where 
the soft mode condenses at the zone center, i.e., qo = 0, 
as in BaTiOa. 

It is well known that the pseudospin and phonon de- 
grees of freedom can be decoupled by introducing dis- 
placed phonon coordinates [3, Efl 

1 Z:^'- 

Qqp = Qqp 7^ E , (10) 

VN 

where 

Yq-r^Y. exp(-z<f • R,) . (11) 

i 

In this so-called adiabatic approach it is implied that the 
time scale for the Ti ion motion is much longer than the 
period of oscillation for phonon modes, so that phonon 
coordinates Q^p adapt adiabatically to any change of the 
pseudospin coordinates Y^r- Thus we obtain the adia- 
batic Hamiltonian 

"^ad — '^dir ^" ^ ^ ^ (^^qpQqpQ — qp ^ ^QP^^QpJ 
qp 

^^EE^r>^-V'- (12) 

ij rr' 

In the last term, the Ti-Ti coupling constant is given by 

qp '^9P 

On a mesoscopic scale, we are not interested in the mi- 
croscopic mechanisms leading to the pseudospin-phonon 
coupling and merely adopt the established functional 
form (lOj . In general, contains the contributions of 
both short- and long-range interactions between the Ti 
and all the other ions. Short-range interactions are ex- 
pected to be primarily responsible for the local potential 
of the Ti ion. Indeed, first-principles calculations indi- 
cate that fcrroelectricity in BaTiO.-j appears as the result 
of hybridization of the Ti-0 bond Therefore, the 

contribution of short-range forces to and hence to 

Kfj" is expected to be dominant, although long-range 

electrostatic dipole-dipole interactions contained in ifj" 
are also needed to establish ferroelectric order 6]. 

The i j part of the last term in Eq. H12(l has the same 
structure as the direct interaction ©. Therefore, Kf^ 



3 



for i ^ j represents an additional interaction between the 
Ti ions at two different sites, which together with the di- 
rect coupling iJJ^ can lead to an order-disorder transition 
into a Tiu polarized state of the Ti subsystem. In gen- 
eral, this interaction involves two different symmetries F 
and F', and is real after being symmetrized with respect 
to the exchange F <-^- F'. Thus we can combine the two 
coupling constants if^^ and KjJ^ into a single coupling 
parameter J^f = if^^ + Kf^ . 

The i ^ j terms in Eq. however, represent a 

constant shift of the local energy and do not contribute to 
the ordering of Ti ions ^3 ■ It can be shown by symmetry 
arguments that Kf^ is zero unless F' = F. 

In the following we will focus on the symmetric part 
of the interaction i^y j = Kf^^. Its Fourier transform is 
given by 



ature. Thus we find 



P 



-E 

q'p 



q'p 



(14) 



The second term is often neglected, however, as shown 
in Ref. [l^ its presence is crucial in order to ensure 
a zero average value of K^. In the q Q limit, we 
have i^or > if X]p l/^pP/'^lp ^ maximum at the 
zone center. For symmetry F = Ti„ this then favors a 
ferroelectric ordering of the Ti subsystem provided that 
Jqt = hr + > 0, implying that (1oti«,s) 0, 
(s = 1, 2, or 3). If, however, the maximum occurs at the 
zone boundary, we can have i^or < 0. For Jgr < 0, the 
ordering is antiferroelectric. It should be stressed that 
in view of relation ij^ the interaction is, in general, 
temperature dependent. 

The off-diagonal coupling J^^ with T' ^ T leads to 
anisotropic interactions, which are assumed to be weaker 
than the isotropic part and can thus be treated by per- 
turbation theory. Below, we will discuss the possibility 
that anisotropic interactions give rise to time-dependent 
random fields acting on the pseudospin variables. 

In a ferroelectric system like BaTiOa, the main contri- 
bution to the coupling K^r will come from the soft TO 
mode with phonon coordinate Qos and frequency wos as 
given by Eq. (j^ . The Ti subsystem will undergo a phase 
transition into an ordered state with nonzero value of the 
pseudospin thermal average (ios) 7^ 0. The transition 
temperature is determined by the relation 



f2 

JOs 



a{Tc - To) 



(15) 



where we have substituted expression for ujq^ and 
written L,o ee {l/N)J2g, l/f'slV'^^-,,- To evaluate L^o, 
we should know the details of the phonon spectrum as 
well as the g-dependence of the coupling. By adding a 
g^-term to wqs in Eq. @ , it can be shown that the leading 
contribution to L^o is a constant independent of temper- 



To 



Mo/k + ^(To-Mo/fc)2+4/2//ca 



(16) 

where Mq = Iq — Lo and we have dropped the indices 
s. Expression H16|l differs from the result of GY by the 
presence of the Mp-term, which is, in general, different 
from zero. We will assume that /q > Lo and < Afo < 
To, implying that T,, > To, i.e., the ordering takes place 
above the stability limit of the unperturbed soft mode. 
It can easily be seen that in the above case, the effect of 
Mo is to shift Tc towards higher temperatures. 



III. SOFT MODE DYNAMICS 

We now consider the response of the soft TO phonon 
mode Q^s to a time-dependent electric field E_^s associ- 
ated with an optic wave in a light scattering experiment. 
Close to the transition, the time scales for the soft mode 
and for the relaxational motion of the Ti ions become 
comparable, and the dynamics based on the adiabatic 
Hamiltonian H12|l is not applicable. The soft mode dy- 
namics is governed by the corresponding part of the orig- 
inal Hamiltonian, namely. 



2 q s q' s —a 



-q's 



2 ^q' \ q's^q's^-q's q's -q' t 



-VW^e fq'sQq-'sY-q-'s - f^E^qsQq, exp(z^i) . (17) 

Here fi = e* j \Jm* is a coupling parameter involving the 
effective charge e* and reduced mass to*. For simplicity, 
we do not include the direct coupling between the light 
vector E^^g and the dipole moment associated with the 
Ti-0 bond. We could, in principle, assume a nonzero 
value of the off-center Ti dipole moment, however, as 
shown by GY its contribution to the Curie constant is 
negligible in other perovskite ferroelectrics and can thus 
be ignored. 

The time evolution of the soft-mode operators Q^s 
and P^s is governed by the Heisenberg equations of mo- 
tion dQ^s/dt = —{i/fi)[Q^s,'Hsiyi\- In contrast, the time 
dependent thermodynamic fluctuations of the variables 
Yis{t) will be assumed to exhibit a pure relaxational 
motion with a single characteristic relaxation time r, 
i.e., we will ignore the possibility of coherent dipole 
moment flips. The corresponding equation of motion 
can be obtained from the classical Langevin model [l5l| . 
which is based on the continuous or 'soft' spin variables 
— oo < Yis{t) < +00 with effective Hamiltonian 

For u — —r ^ 00 one recovers the discrete limit = 1 
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(cf. Eq. (O). The Langevin equation of motion is 



dt 



(19) 



where the Langevin noise ^is{t) is a Gaussian random 
variable with zero mean and variance 



(20) 



Introducing the Fourier components Y^s{uj) etc., we ob- 
tain the Hnearized equations of motion 



(21a) 
(21b) 

-(r - f3I^, - ^^s)Y^s + PVNf^.Q^, + ^^,{lo) . 

(21c) 



t^Qqs ~ Pqs •) 

2 1 

V Jy 



Here S^s is the self energy which can, in principle, be cal- 
culated by a diagrammatic expansion involving the pa- 
rameter u and the pseudospin-phonon coupling /g-^.. In 
the following we will ignore the frequency dependence of 
Tt^s in the soft-mode regime. 

We can now introduce the static pseudospin response 
Xqs = {SY^s/S^qs), i-e.. 



Xqs 



(22) 



and redefine the relaxation time by writing t^s = rxqs/ P- 
The solution of the above equations can then be ex- 
pressed in the form 



and 



V I \ fqsXqs / \ 

Yqs{uj) = ^ , . Qqs{uJ) , 



nEq-, (23) 



(24) 



where the random force term has been averaged out. 
Eq. H23(l is similar, but not identical to the expression 
given by GY, which contains an extra term ituTqs in 
the denominator, describing the damping of the TO soft 
mode due to lattice anharmonicity. Here we absorb this 
damping into the ilut^s term by a proper redefinition the 
relaxation time t^s- 

We now consider the q = case corresponding to a 
TO soft mode at the zone center. Dropping the sub- 
scripts s and introducing the far-infrared dielectric re- 
sponse xq{^) = Qo{^^)/{f^EQ) we get 



xq(^) 



1 + iujTo 



(25) 



The contribution to the corresponding dielectric function 
is given by 



£0^0 



■Xq(w) 



(26) 



where vq is the unit cell volume. 

In the absence of pseudospin-phonon coupling, /q — > 0, 
the static response xq(0) diverges at Tc = To in view 
of Eq. ©. For /o ^ 0, the critical temperature is ob- 



tained from the equation ujq - 
Eq. from the relation 

a{T - To 



'/oXo = 0, or after applying 



(r - (3Io - So) 



(27) 



Returning to the discrete limit Y^^ = 1 we notice that the 
last result will be equivalent to Eq. (|15|l provided that 
r — Eq ~ 1 -|- fiLo- Thus, in the above limit the critical 
temperature Tc derived from the soft-mode dynamics will 
be precisely equal to the static value (|15|l . This means 
that the hybrid soft mode becomes unstable at the static 
ordering temperature of the Ti subsystem. It should 
be added that the technical reason for using the soft spin 
variables and Langevin dynamics is the fact that only in 
this formalism we were able to properly take into account 
the presence of the second term in Eq. (|14|l . 

To illustrate the frequency dependence of Xq(^) in- 
troduce dimensionless parameters /o —> foTo/VkTo and 
a —> qToTq, and choose their representative values as 
fo — 0.4 and a ~ 0.145, respectively. Also, we fix the 
value of the parameter Mq at Mo/kTo — 0.5 and rescale 
the frequencies as — > luto- In all numerical calculations 
we furthermore set tq = const., i.e., we ignore the non- 
critical temperature dependence of the relaxation time. 
We can then evaluate the spectral function 



UJ 



(28) 



which is related to the Raman scattering intensity. Here 
the imaginary part of the response (|25l) is defined by 
Xq{^) = XQ(w)-iXQ(w). In Fig. 1, we plot /q(w) at five 
different temperatures T/T^. Far from T^,, the spectrum 
is characterized by a single peak, the position of which 
moves towards zero frequency as T — > Tc in accordance 
with the expected soft-mode behavior. 

The hybrid soft-mode frequency approximately corre- 
sponds to the position of the peak in Iq(lu). To obtain 
an unambiguous definition it is necessary to consider the 
complex poles of the response function (|22ll . After intro- 
ducing a new variable z = ituTo, the denominator takes 
the standard form of a cubic polynomial with real coeffi- 
cients 



3t 2, 22 i/2 j-2 \ 2 

+ z +ujgT z + (cjg - /oXo)to 

= (Z - Zl){z - Z2){Z - Z3) = . 



(29) 



The complex zeros 2:„, n — 1, 2, 3, are given by the Car- 
dano formula. The corresponding analytic expressions 
for z„ are cumbersome, but can easily be evaluated nu- 
merically. They then lead to three complex frequencies 
— Xn + iyn having the following properties: (i) One 
of the frequencies is always imaginary, say, uji = iyi; (ii) 
The remaining two solutions are, in general, complex and 



5 



such that always = X2] (iii) If X2 = 0, then ^ y2, 
however, for X2 ^ we have y-^ = y2', (iv) In general, 

J^nUn = l/'^o, and (v) ujiLl>2 +^^2^3 + uj^^i = -^o- 
The response H25|l thus becomes 



UJ — IVq 



(30) 



where vq = 1 /tq and the resonance frequency flo is given 
by 



"0 



-UJ2i^3 =X2+ 2/22/3 , 



while the damping parameter is 

7 = 2/2 + 2/3- 



(31) 



(32) 



The term i{x2yz + 2:32/2) formally appearing in the prod- 
uct W2W3 is zero in view of the properties (ii) and (iii) 
above. It further follows from (iv) that yi — i^q — 7, 
whereas from (v) we obtain a relation between Qq and 7, 
namely, 



■0 + 7(^0 - 7) = ' 



We can now rewrite the spectral function 
form 



A 



B-Aoj^ 



Vl 



172)2 +CJ272 



(33) 



in the 



(34) 



where A = + yf- 2/17) and B = {2^1 + yf-T 

The first term represents a central peak associated with 
the relaxational motion of the Ti ions, whereas the second 
term describes a resonance occurring at Qq. It should 
be noted, however, that because of the —Aoj'^ term in 
the numerator the last expression deviates from the stan- 
dard spectral function for the damped harmonic oscilla- 
tor. Furthermore, the —Aoj'^ term partially compensates 
the contribution of the relaxational term in Eq. (|34|l . so 
that a separate central peak cannot be discerned at tem- 
peratures of interest, as it is evident from Fig. 1. 

In Fig. 2, we plot JIq as a function of temperature 
with the same parameter values as in Fig. 1. For /o 
smaller than a limiting value /™, where = 0.6721 for 
the present choice of representative parameter values, the 
temperature dependence of fig near Tc is quasi linear. 



"2)A. 



nl~a'iT-Tc) 



(35) 



where a' = a'(/o) with a' > a and Tc > Tq. Thus we 
are dealing with a hybrid soft mode, which condenses at 
the instability temperature Tc given by the static critical 
temperature IjlSfl . 

The temperature dependence of the damping constant 
7 is shown in Fig. 3. For /o < /™ we have < 7 < I'o, 
whereas for fo > we can see that 7 ^ i^o as T ^ Tc. 
It can be shown that in the hybrid soft mode regime fo < 
/™ the reduced damping constant 70 = j/i^o diverges on 
approaching Tc from above as 70 (T — Tc)~^/2_ The 
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FIG. 1: Spectral density of far-infrared dielectric re- 
sponse, plotted on a logarithmic vertical scale, calculated for 
pseudospin-phonon coupling strength fo = 0.4 and for several 
values of temperature, as indicated. 



temperature dependence of 70 is shown in the inset of 
Fig. 3. 

The static susceptibility as obtained from Eq. (|3()|l is 
given by 



xq(0) 



1 



rof7g(i^o - 7) 



(36) 



and is found to diverge at Tc as Xq(0) ~ (T—Tc)^^ for all 
values of /q. As already stated, for fo > /q" the damp- 
ing parameter 7 approaches the value 7 — > I'o as T — > Tc, 
and thus Xq(0) diverges even for flo{Tc) ^ 0. The tem- 
perature dependence of the inverse static response Xq 
is displayed in Fig. 2. For T < we formally have 
ilo < and the system becomes unstable. In principle, 
it should be stabilized by including anharmonic terms 
into the phonon Hamiltonian © as well as a coupling to 
elastic strains. 

For fo > /™, ilo is always nonzero and strictly speak- 
ing the hybrid soft-mode concept is not applicable in 
that case. The corresponding limiting value of the in- 
stability temperature corresponding to the representative 
parameter values used in numerical calculations here is 
T^ /To = 2.533. This is a fairly large value, suggest- 
ing that any real system is likely to be in the soft-mode 
regime < /o < /^'. 
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12 3 4 



FIG. 2: Resonance frequency Qo as a function of temperature 
for several values of coupling strength fo (solid lines). For 
fo < fir ~ 0.6721 a soft-mode behavior is observed. Also 
shown is the temperature dependence of the inverse static 
response (dashed lines). 



1.2 




12 3 4 



FIG. 3: Damping parameter 7 as a function of temperature 
for several values of coupling strength fo. Inset: Reduced 
damping parameter 70 — ^/Q,o for the same values of fo. 



IV. DISCUSSION 

To determine the type of the phase trasition occur- 
ring at T — Tc, we introduce the order parameter as the 
thermal average of the Ti pseudospin variable P = (Yq), 
which is related to (Qq) through Eq. at w = 0. We 
can then write down a Landau-type free energy 

HP) - I Xy{0)-'p' + i + i + . . . , (37) 

where the static pseudospin response xy{0) is related to 
the hybrid soft-mode response H25|l and the rigid static 
pseudospin response H22() . namely, 

Xy{0) = XQ{0)^^^^ (7^) • (38) 
^0 - JoXo V io Xo / 

As noted in deriving Eq. (|27|1 . the first denominator van- 
ishes as T Tc, while in the same limit the last factor 
approaches the value (^q/JqXo 1- Eq. H38|l is reminis- 
cent of the Lyddane-Sachs- Teller relation for the static 
dielectric response. It follows that the coefficient of the 
quadratic term in !F{P) vanishes at T,. as ^ — Tc). 
The coefficients h and c can easily be determined within a 
molecular field approximation, however, it turns out that 
this approximation would predict a second order phase 
transition, in disagreement with observations in BaTiOa. 
To derive the correct value of the coefficient b in the free 
energy we would have to include, for example, a coupling 
between the soft mode and elastic strain or equiva- 
lently, between the strain and the Ti pseudospin variable. 
The appearance of a nonzero value of the order parameter 
will thus be accompanied by a macroscopic deformation 
of the lattice with tetragonal symmetry. In practice, h 
and c are often considered as p henomenological param- 
eters. Hatta and Ikushima |l7l | applied an expression of 
the above form to analyze the measured heat capacity of 
BaTiOa at constant electric field. They used the value 
Tc = 383 K (in our notation) and found that 6 < and 
c > 0, implying a first order phase transition which takes 
place at Tc — 398 K (or at ~ 408 K, depending on the 
sample preparation method). They also determined the 
jump in the heat capacity at Tc of the order AC ~ 0.19 
k, from which it was concluded that the Ti ion lies in 
a single minimum potential, in apparent disagreement 
with the off-center Ti picture. According to the present 
model the number of equilibrium positions of Ti is greater 
than one, but their actual number is only relevant for 
the definition of the Ti pseudospin variables. If we re- 
strict the discussion to the adiabatic Hamiltonian H12|l . 
the phase transition appears to be of order-disorder type. 
However, the above approach involving hybrid soft-mode 
dynamics — which leads to the same critical temperature 
as the static approach — clearly has the characteristics of 
a displacive transition, but with a simultaneous order- 
ing of the Ti subsystem as an additional order-disorder 
feature. 

Evidence of a TO soft mode with Ti„ symmetry in 
BaTiOa has been found by Vogt et al. [3| using the hyper- 
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Raman scattering technique. The dielectric function has 
been fitted by a classical single oscillator dispersion for- 
mula, which is equivalent to the second term of Eq. H34|l 
without the cj^-term in the numerator. One of the im- 
portant results of Ref. 2] is that the frequency of the 
soft mode tends to zero close to the transition tempera- 
ture and no saturation of the frequency at a finite value 
reported earlier occurs ,181] , in agreement with the predic- 
tions of the present model. Close to Tc, the temperature 
dependence of the soft-mode frequency can be described 
by the relation ^ with a' ~ 1.4 x lO^^ s-^K"!. The 
relative damping constant 7/O0 was found to exceed the 
value of 2 at all temperatures, whereas in the present 
model this is true only close to Tc. 

It has been shown by NMR methods 5] that in the high 
temperature phase of BaTiOs the unit cells are tetrago- 
nally distorted, although the overall macroscorac symme- 
try is cubic. This agrees with Wada et al. who de- 
termined the symmetry by Raman scattering and found 
that the microscopic symmetry was P4mm both above 
and below Tc. Since the orientation of the tetragonal 
axis varies across the crystal, the macroscopic symmetry 
is Pm3m above the transition and P4mm below Tc- This 
can be accounted for by the present model if we return 
to the Ti-Ti coupling terms (0) and H13() . and allow for 
off-diagonal interactions with F' ^ F. This then gives 
rise to extra terms the equations of motion H21c|l for Y^g , 
e-g-, 

iutY^ ^ -{r- 131^ - S,-r)iV + PVNf^^Q^ 

+PEr'^rJr'Y?r'+^?ri^) ■ (39) 

On the time scale of the TO soft mode and of the Ti 
relaxational mode — Y^g the variables F^t' are much 
slower and appear to be "frozen" in a given configura- 
tion. Thus the sum over F' will play the role of a random 
variable = X^r'^^r -^f^ ^qr'{t) analogous to the ran- 
dom electric field in dipolar glasses. Instead of spatial 
randomness, however, we are dealing here with a tempo- 
ral disorder, which appears to be "quenched" on the time 
scale of YqT (t) ■ This field will give rise to a slowly varying 
deformation of Ti^^s symmetry, which is experimentally 
observable in both NMR and Raman experiments. Since 
the orientation of the deformation axis s varies both in 
space and time, the average symmetry of the system re- 
mains cubic. 



V. CONCLUSIONS 

Recent experiments provided new evidence that the 
Ti ion in the cubic phase of BaTiOa ^ and other ox- 
ide ferroelectrics |2^ occupies one of the eight off-center 
positions in the unit cell along the [111] directions. We 
have presented a simple mesoscopic model of BaTiOs, as- 
suming that the symmetrized occupational probabilities 
for the Ti sites or pseudospins are linearly coupled with 
the normal coordinates of lattice vibrations. On the time 



scale of Ti intersite jumps, lattice modes act as fast vari- 
ables, which can adapt instantaneously to any change of 
the pseudospin configuration. In this adiabatic approx- 
imation, the pseudospin-phonon coupling gives rise to a 
static phonon mediated effective Ti-Ti interaction. The 
Ti subsystem thus undergoes an order-disorder transition 
into a polarized state of Ti„ symmetry, and the critical 
temperature Tc is determined by a sum of the contribu- 
tions from the direct Ti-Ti coupling and the pseudospin- 
lattice part. The leading contribution is due to the TO 
soft mode, the frequency of which tends to zero at the 
instability temperature Tq. In general, Tc > Tq, where 
the shift of Tc depends both on the direct interaction as 
well as the pseudospin-phonon coupling. 

Close to the stability limit for he TO soft mode, the 
time scales for the Ti and lattice motion become compa- 
rable and a dynamic treatment becomes necessary. By 
combining the classical Langevin equations of motion for 
the relaxation of the Ti pseudospins and the quantum 
Heisenberg equations for the TO soft mode coordinate, 
we have obtained the response of the system to an oscil- 
lating electric field, typically in an optical scattering ex- 
periment. The spectral function of this response consists 
of a central peak and a resonance, belonging to a new hy- 
brid soft mode with Tiu symmetry, which describes the 
in-phase motion of the TO soft optic mode and the re- 
laxational decay of the symmetrized Ti pseudospin vari- 
able. The frequency of the hybrid soft mode tends to 
zero at the static critical temperature Tc. The central 
peak component is determined by the relaxational mo- 
tion of the Ti subsystem and can be regarded as a typical 
order-disorder feature of the total response, in contrast to 
the resonance which shows a mixed displacive and order- 
disorder character. 

The TO soft mode in BaTiOa has been observed exper- 
imentally by hyper- Raman scattering 0, and has been 
described by a damped harmonic oscillator, similar to 
our resonance term in the dynamic response. The soft 
mode frequency was found to extrapolate to zero close to 
the thermodynamic transition temperature Tc- 

The frequency of the hybrid soft mode depends on the 
Ti-phonon coupling strength and exists only in a finite 
interval < /o < /™, where the value of /™ depends 
on the model parameters. Outside this interval, the soft 
mode concept is not applicable, however, the static re- 
sponse always diverges at Tc. 

To calculate the free energy of the system we would 
have to include a coupling between the soft mode and 
lattice strains into the model. This can be done in a 
phenomenological approach "I6] , which describes the first 
order phase transition from the high temperature phase 
with macroscopic cubic symmetry to a tetragonal low 
temperature phase. The corresponding Curie tempera- 
ture Tq is generally higher than the instability temper- 
ature Tc associated with the condensation of the hybrid 
soft mode. 
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